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ABSTRACT 

We use a sample of 7 starburst galaxies at intermediate redshifts {z ~ 0.4 and z ~ 0.8) with observa- 
tions ranging from the observed ultraviolet to 1.4 GHz, to compare the star formation rate (SFR) estima- 
tors which are used in the different wavelength regimes. We find that extinction corrected Ha underesti- 
mates the SFR, and the degree of this underestimation increases with the infrared luminosity of the galax- 
ies. Galaxies with very different levels of dust extinction as measured with SFRir/SFR(Hq;, uncorrected 
for extinction) present a similar attenuation A [Ha], as if the Balmer lines probed a different region of 
the galaxy than the one responsible for the bulk of the IR luminosity for large SFRs. In addition, SFR 
estimates derived from [O ii]A3727 match very well those inferred from Ha after applying the metallicity 
correction derived from local galaxies. SFRs estimated from the UV luminosities show a dichotomic 
behavior, similar to that previously reported by other authors in galaxies at z < 0.4. Here we extend 
this result up to z 0.8. Finally, one of the studied objects is a luminous compact galaxy (LCG) that 
may be suffering similar dust-enshrouded star formation episodes. These results highlight the relevance 
of quantifying the actual Lir of LCGs, as well as that of a much larger and generic sample of luminous 
infrared galaxies, which will be possible after the launch of SIRTF. 

Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: starburst — galaxies: stellar 
content 



1. INTRODUCTION 

The observational efforts devoted during the last years 
to measure the cosmic star formation history is clearly 
manifest in the growing number of works in the field. As 
a result of the advent of the Hubble Space Telescope and 
of the 8-10 meter class telescopes, the rest-frame ultravio- 
let (UV) star formation history per unit comoving volume 
initially derived up to z ^1 (Lilly et al. 1996) was ex- 
tended up to z - 4 (Madau et al. 1996; Steidel et al. 1999). 
Naturally, all these works rely on the use of the available 
star formation rate (SFR) diagnostics: nebular emission 
lines Hke Ha (Kennicutt 1992a, 1998; Gallego et al. 1995; 
Tresse & Maddox 1998; Glazcbrook et al. 1999; Moorwood 
et al. 2000; Hopkins et al. 2000) and [O ii] (Hogg et al. 
1998; Aragon-Salamanca et al. 2002); UV continuum lu- 
minosities (Lilly et al. 1996; Madau et al. 1996; Treyer 
et al. 1998; Cowie, Songalia & Barger 1999; Steidel et al. 
1999; Sullivan et al. 2000, 2001); far-infrared luminosities 
(Rowan-Robinson et al. 1997; Blain et al. 1999; Flores et 



al. 1999; Chary & Elbaz 2002); and 1.4 GHz radio lumi- 
nosities (Condon 1992; Cram et al. 1998). Different cosmic 
star formation histories have been advocated from these 
various studies which can only be reconciled if the effect 
of dust extinction is accounted for. 

In the local universe, galaxies radiate the bulk of their 
luminosity in the stellar regime, i.e. below A ^ 5 /im. On 
average, only about one third of their stellar light is ab- 
sorbed by dust and re-emitted in the dust regime from the 
mid infrared (MIR) to the sub-millimeter (Soifer & Neuge- 
bauer 1991). In the distant universe, the reverse seems to 
take place as suggested by the combined results of the 
COsmic Background Experiment (COBE), the ISOCAM 
and ISOPHOT instruments on board the Infrared Space 
Observatory (ISO) and the sub-millimeter camera SCUBA 
at the JCMT. Number counts at 15 /im (Elbaz et al. 1999), 
at 90 fim and 170 /xm (Efstathiou et al. 2000; Kawara et 
al. 2000; Matsuhara et al. 2000; Dole et al. 2001), and at 
450 and 850 /im (Small et al. 2002, and references therein) 
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revealed the presence of an excess of faint galaxies by one 
order of magnitude in comparison with expectations based 
on the local properties of galaxies. Such excess can only 
be explained if galaxies emitted a larger fraction of their 
luminosity in the dust regime in the past. This conclusion 
is reinforced by the detection of a strong cosmic infrared 
background (CIRB) above A ~ 100 /xm (Puget et al. 1996; 
Hauser et al. 1998; Fixsen et al. 1998; Lagache et al. 1999, 
2000) the bulk of which appears to be produced by star- 
bursting galaxies radiating the majority of their light in 
the IR above A ~ 5 /xm (Elbaz et al. 2002; Chary & Elbaz 
2002). 

Two major questions remain unanswered: 

- How consistent are the different SFR indicators com- 
monly used in the literature? 

- Do local and distant galaxies exhibit a similar be- 
haviour? 

Studies of local galaxies have shown that even after cor- 
recting for dust extinction, a significant discrepancy re- 
mains between SFRs derived from UV and Ha measure- 
ments compared with those obtained from FIR and radio 
luminosities, and that this discrepancy is worse for galax- 
ies with higher SFRs (Hopkins et al. 2001; Poggianti & 
Wu 2000). Hopkins et al. (2001) and Sullivan et al. (2001) 
report an increase of the F(Ha) /F(H/3) ratio, used to mea- 
sure extinction, with the SFR of galaxies. They suggest 
to use the fit to this correlation in order to account for 
dust extinction when H/3 is not available. However, Buat 
et al. (2002) indicate that the apparent link between ex- 
tinction and star formation is likely not real, but only the 
result of a dispersed correlation between dust extinction 
and luminosity. 

In a recent paper, Rigopoulou et al. (2000) studied a 
sample of 12 galaxies with 0.4 < z < 1.4, detected with 
ISOCAM in the Hubble Deep Field South (HDF-S). They 
showed that the discrepancy found between optical-UV 
and FIR-radio SFRs in the local universe is also present 
at high z. However, the Balmcr decrement was not avail- 
able for these galaxies and the authors discussed the effect 
of dust extinction assuming an average correction factor 
of 4 for SFR(Ha). 

In the present paper, we propose to address the two 
previous questions using a test sample of seven galaxies 
at 2: ~ 0.4 and 0.8, five of which have MIR flux densities 
from which we derived FIR luminosities assuming that the 
correlation betwecni AIIR and FIR luminosities observed in 
the local universe remains valid at these redshifts (Elbaz et 
al. 2002; Chary & Elbaz 2002). The galaxies were selected 
from their MIR flux density in order to sample SFRs in 
the range from 2 to 160 Mq yr^^ and morphologies ranging 
from a normal spiral to merging galaxies. 

In this work we extend to higher redshifts the compari- 
son between different SFR estimators that have been used 
for samples of local galaxies (e.g., Sullivan et al. 2001). For 
that purpose, we compare classical SFR indicators such 
as the nebular emission lines [O 11] and Ha, or the UV 
continuum, that we correct for dust extinction using the 
Balmer decrement, with the FIR and radio luminosities 
which are unaffected by extinction. This has been made 
possible with the comissioning of high resolution optical 
and near-infrared (NIR) spectroscopy on the Keck tele- 
scopes. High resolution spectroscopy is required not only 



to avoid the contamination from atmospheric OH emis- 
sion lines but also to resolve the [N ll]AA6549, 6583 A and 
Ha emission lines. After the presentation of the observa- 
tional data in Sect. 2, the fluxes in these three lines will 
be combined in Sect. 3 with the [O ill] and H/3 lines to lo- 
cate galaxies in diagnostic diagrams that allow separating 
between galaxies whose optical luminosity is mainly due 
to star formation, from galaxies where an Active Galactic 
Nucleus is the dominant source of emission. 

In Sect. 4 we obtain imiltiwavclcngth SFR estimates for 
the galaxy sample using the available flux data. Addi- 
tional galaxy parameters are examined in Sect. 5, and a 
summary of the main properties of each object is presented 
in Sect. 6. In Sect. 7 we compare and discuss the results 
obtained using the different SFR estimators. Finally, we 
present the conclusions in Sect. 8. 

Unless otherwise indicated, throughout this paper we 
will assume Hq = 70 km Mpc~^, Qm = 0.3, and 
= 0.7. 

2. THE DATA 

2.1. Galaxy Sample and Observations 

The sample of dusty starburst galaxies was chosen ac- 
cording to the availability of ISOCAM-15 /xm luminosities 
(or, at least, upper limits) for each object. The galaxies 
were also constrained to have a spectroscopic redshift ac- 
curate enough to insure that the Ha line would be free of 
contamination due to night sky lines. With those criteria 
in mind, we have selected 4 galaxies with a mean red- 
shift {z) ~ 0.46 in the Hubble Deep Field North, HDF-N 
(Williams et al. 1996), and 3 galaxies with (z) ~ 0.80 
in the Groth Strip Survey, GSS (Groth et al. 1994; Koo 
et al. 1996). The whole set of galaxies is listed in Ta- 
ble 1 and displayed in Fig. 1. Interestingly, the two galax- 
ies undetected by ISO have a compact morphology (hd2- 
264.2 and GSS084.4515). The ISO sources are face-on spi- 
rals (hd4-656.1, hd4- 795.111, and GSS084_4521), a disky 
galaxy (hd2-264.1), and a colliding system (GSS073_1810) 
which looks like a higher redshift analogue of the local 
Antenae system (NGC 4038 and NGC 4039). In order to 
facilitate the recognition of the individual objects in some 
of the flgures of this paper, we have labeled the symbols 
with the letters S, D, C and A, corresponding to their mor- 
phologies, i.e. Spiral, Disky, Compact and Antenae-like, 
respectively. 

According to their IR luminosities (computed as de- 
scribed in Sect. 4.3 and listed in Table 9), the HDF-N 
galaxies are not luminous enough to be classified as lumi- 
nous infrared galaxies (LIRG, iiR, > lO^^L©), whereas the 
contrary is true for the GSS objects. Note that although 
GSS073_1810 could be marginally considered as an ultra- 
luminous infrared galaxy (ULIRG; Lir > lO^^L©), this 
object is constituted by two colliding galaxies (see Fig. 1). 
A more detailed characterization of the galaxy sample is 
presented in Sects. 5 and 6. 

For all galaxies in our sample, we have obtained opti- 
cal and near-infrared spectroscopy. The observations were 
carried out with the W. M. Keck Observatory telescopes 
and 3 different instruments over 5 observing runs (see Ta- 
ble 1). For the relatively low redshift HDF-N galaxies, 
we employed the Echelle Spectrograph and Imager, ESI 
(Sutin 1997; Epps & Miller 1998; Sheinis et al. 2002), 
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which covers the full wavelength range from 3900 A to 
11000 A in one configuration, and with a spectral res- 
olution FWHM ~ 0.8 A at observed A(Ha) ~ 9600 A. 
However, for the more distant GSS galaxies, initially only 
observed in the spectral range 4500 9000 A with the Low- 
Resolution Imaging Spectrometer, LRIS (Oke et al. 1995) 
(FWHM ~ 3 A), we expanded the spectral window by 
using the near-infrared echelle spectrograph NIRSPEC 
(McLean et al. 1998) in the J band (N2 and N3 filters, 
total spectral range 1.089-1.375 /xm, FWHM - 8 A). In 
general the observations were carried out with airmasses 
ranging from 1.09 to 1.56. 

2.2. Data Reduction and Error Handling 

The reduction of the LRIS data (runs 1, 2 and 3) is as de- 
scribed in Simard et al. (1999). The absolute flux calibra- 
tion used observations of the spectrophotometric standard 
BD+33 2642 (runs 1 and 3) as the primary calibrator and 
the spectra of two stars, GSS063_0445 and GSS093_2453 
(runs 1, 2 and 3), used to align the image masks, as sec- 
ondary calibrators. We find that the response curve of 
LRIS is very stable on a run-to-run basis, with small fluc- 
tuations < 5%, which translate into negligible variations 
in the star formation rates once other sources of error and 
uncertainties are considered. 

The ESI (run 4) and NIRSPEC (run 5) spectra have 
been reduced using a suite of our own programs based on 
the REDUCEME package (Cardiel 1999), and a new FOR- 
TRAN package, XNIRSPEC, specially developed for the re- 
duction of the near-IR data. Error frames have been com- 
puted directly from the raw images, using the read-out 
noise and the gain of each detector. The reduction in 
parallel of these error frames together with the data im- 
ages guarantees the correct propagation of errors due to 
arithmetic manipulation throughout the reduction proce- 
dure (Cardiel et al. 1998, 2002). In this sense, each fully 
processed spectrum is accompanied by an associated error 
spectrum which contains the random error in each pixel 
due to photon statistics and read-out noise. 

One of the major problems encountered when measuring 
spectral features in the near-IR is the proper removal of 
the OH emission lines (Rousselot et al. 2000) and the com- 
pensation for telluric absorption due to water vapour and 
other molecules (Stevenson 1994; Chmielewski 2000). Sys- 
tematic variations in the spectral direction (e.g. sampling 
aliasing, inacurate geometric distortion corrections), and 
in flux calibration (e.g. even small flatfleld biases, residual 
fringing, or slit width inhomogeneities), pose a difficult 
challenge to the data reduction. Trying to minimize part 
of these problems, the NIRSPEC observations were car- 
ried out nodding the telescope, changing the position of 
the galaxies on the slit, in order to feed the data reduction 
pipeline with the subtraction of consecutive images. Even 
when following this procedure, the unpredictable residu- 
als of each sky line (either in intensity or sign) had to be 
subtracted like normal sky lines in a subsequent step. 

Since the spectra obtained in the detector plane of ESI 
and NIRSPEC are distorted, the image rectification prior 
to the removal of sky lines generates aliasing artifacts (the 
undersampled sky line residuals cannot be properly rec- 
tified), and, in addition, introduces a correlation between 
adjacent pixels in the error frames. In order to alleviate 



both effects, we mapped the image distortions in the spa- 
tial and spectral directions by using bivariate polynomial 
transformations (Wolberg 1990). These mapping func- 
tions were determined by fitting the spectrum of a bright 
star, observed at diff'erent positions along the slit, and to 
the sky lines. 

For the ESI observations, the removal of the telluric ab- 
sorption, as well as the absolute fiux calibration, used 
the high signal-to-noise ratio (S/N per A > 150) al- 
most featureless spectra of the spectrophotometric stars 
(BD-l-28 4211, Wolf 1346 and Feige 34). In the NIR- 
SPEC run only one spectrophotometric star was observed 
(BD-(-28 4211), and the removal of telluric features was 
performed using the observation of V986 Oph, a hot 
(BOIII) and bright (V = 6.15 mag) star. 

All the fiux calibrated spectra are displayed in Figs. 2 
and 3. 

2.3. Measure of emission line fluxes 

The spectra have been corrected for atmospheric extinc- 
tion using the Mauna Kea values available on the WWW 
page of the UKIRT (United Kingdom Infrared Telescope). 
The small corrections for galactic extinction used the dust 
maps of Schlegel et al. (1998) and the galactic extinction 
curve of Fitzpatrick (1999). 

With the aim of correcting the measured Balmer emis- 
sion line fluxes for the expected underlying stellar absorp- 
tion, we have applied the criterion followed by Calzetti et 
al. (1994), who, based on the extensive analysis of giant ex- 
tragalactic HII regions by McCall et al. (1985), employed 
EWabs(HQ!) = EWabs(H/3) = 2 A for galaxies with unde- 
tected or uncertain H7 flux. 

In order to follow a homogeneous procedure in the de- 
termination of line luminosities, line fluxes have been de- 
rived from Gaussian fits to the corresponding emission 
lines. It is worth noting that these fits not only match the 
observed spectra within the error bars available in each 
pixel, but they also help minimizing the uncertainties due 
to the choice of integration limits on both sides of the 
emission lines. Flux errors are also easily computed in this 
way through numerical simulations with bootstrapped er- 
ror spectra. All the measured emission lines are displayed 
in the insets of Figs. 2 and 3 (see caption for details). Their 
values are listed in Table 2. Considering that the Ha line 
has been observed in spectral regions densely populated by 
sky lines, it is not surprising that in some cases large resid- 
uals prevented fitting Gaussians to the entire wavelength 
domain of the [N 11] AA6549, 6583 lines. In these cases (and 
also when no clear signature of these lines was detectable), 
since the expected location of the lines is very well con- 
strained by the centroid of the Ha line, restricted fits were 
performed to the line regions free of contamination (see 
caption of Table 2 for details). Restricted fits have also 
been applied to other emission lines when required. 

Intrinsic reddenings, parameterized through the color 
excess in the nebular gas, E(B— V)gas, were determined us- 
ing the measured Balmer decrement, and assuming typical 
gas conditions, i.e. A^c = 10^ cm~"^, Tc = 10^ K, recom- 
bination case B, and the flux ratio Ha/H/3=2.86 (Oster- 
brock 1989). Since the extinction properties of the Milky 
Way (Savage & Mathis 1979; Seaton 1979; Cardelli et al. 
1989; O'Donnell 1994; Fitzpatrick 1999), the Small Magel- 
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lanic Cloud (Prcvot ct al. 1984; Bouchet et al. 1985), and 
the Large Magellanic Cloud (Howarth 1983; Fitzpatrick 
1985) are similar in the optical and in the near-IR spectral 
ranges, it is common to assume that the selective extinc- 
tion curve from any of these objects can be safely used to 
correct the emission lines for internal reddening (Calzetti 
1997). For this reason, here we adopt the average curve for 
Galactic extinction published by Fitzpatrick (1999). The 
E(B— V)gas values for our galaxy sample are given in the 
last column of Tabic 2. The; quoted errors correspond to 
the propagation of random errors in the Ha and H/3 fluxes. 

Finally, to estimate the impact of the slit widths and 
position angles in the measured emission line fluxes, we 
computed the fraction of light passing through the slit 
for each observation, using the HST images; these frac- 
tions and associated errors, are listed in Table 3. The 
errors were calculated using numerical simulations which 
assume random positioning errors when placing the slits 
over the targets. The simulations were performed after 
isolating each galaxy from its neighbors using SExtractor 
(Bertin and Arnouts 1996), with the configuration param- 
eter CHECKIMAGE_TYPE set to OBJECTS. Although those 
slit-loss estimates must be taken with caution, it seems 
unlikely that we are losing more than 50% of the light in 
any case. 

2.4. Additional Flux Data 

Additional photometric data, shown in Tables 4, 5, 
and 6, were used as consistency checks on the absolute flux 
calibration over the entire wavelength range. These data 
were also used to estimate the importance of the aperture 
correction when measuring the emission line fluxes (see 
Figs. 2 and 3), and to derive SFRs based on ultraviolet, 
infrared and radio luminosities. 

The HST magnitudes for the HDF-N galaxies (Table 4) 
correspond to the fluxes published by Fernandez-Soto et 
al. (1999). 

For the GSS galaxies (Table 5) we list the U, B, R and / 
band magnitudes determined by Brunner et al. (1999) us- 
ing photometric data obtained with the Prime Focus CCD 
(PFCCD) camera on the Mayall 4 m telescope at Kitt Peak 
National Observatory. In addition, we have computed Veoe 
and /814 band magnitudes using SExtractor (Bertin and 
Arnouts 1996) applied on the HST images (Groth et al. 
1994). We have also compared the KPNO magnitudes 
with their analogous HST measurements. In this sense, 
a small offset of 0.15 mag has been detected between the 
KPNO / band magnitudes and the HST Isi4 data (see 
Fig. 4). We do not make a similar correction for the 
KPNO R band magnitudes because the overlap in wave- 
length with the 1^606 bandpass is small and the compari- 
son between both magnitudes has a large scatter. Finally, 
we also derived a direct measurement of the J band mag- 
nitude for the pair of interacting galaxies GSS073_1810. 
In this case, we made use of the images obtained with the 
sUt-viewing camera SCAM of NIRSPEC through the N3 
filter (covering the range 1.143-1.375 microns). Although 
the total effective exposure time only lasted 80 seconds, it 
was enough to get the absolute flux with a relative random 
uncertainty of ~ 25%. 

The MIR flux densities (Table 6), S^{15iJ.m), are the 
ISOCAM measurements with the LW3 filter published by 



Ausscl et al. (1999) and Flores ct al. (1999) for the HDF-N 
and the GSS galaxies, respectively. Non-detections are in- 
dicated by an upper limit, that corresponds to the sen- 
sitivity threshold of the surveys. The flux published by 
Ausscl ct al. (1999) for hd2-264.1 integrates the emission 
of four blended sources. The flux value given here is the 
result for only this galaxy, obtained after deconvolving the 
different objects (H. Aussel 2001, private communication). 
The fluxes for GSS073_1810 and GSS084.4521 are also re- 
vised values (after improving the flux calibration through 
Monte Carlo simulations; H. Flores 2001, private commu- 
nication). 

The radio flux densities (last column in Table 6), 
S'y (radio), arc the 8.5 GHz measurements from Richards 
et al. (1998) for the HDF-N galaxies, and the 5.0 GHz 
data from Fomalont et al. (1991) for the GSS objects. We 
also quote the sensitivity limit of the considered survey for 
the undetected galaxies. Finally, to be consistent with the 
estimation of the MIR flux for hd2-264.1, we consider that 
the radio flux in this case also integrates the combined val- 
ues of the blended sources. Thus, the quoted radio flux is 
the same fraction of the total radio flux. 

3. CONSTRAINING THE AGN CONTRIBUTION 

Before calculating the SFRs from the collected multi- 
wavelength data, we must consider the possible contami- 
nation of our galaxy sample by AGNs. The fact that QSOs 
and ULIRGs exhibit similar total luminosities and space 
densities led Sanders et al. (1988) to suggest that ULIRGs 
may constitute the dust-enshrouded phase of QSO forma- 
tion. In this sense, strong interactions and mergers are the 
mechanisms responsible for the funneling of gas into the 
central regions of galaxies, which subsequently may feed 
both star formation and AGN activity. Although some as- 
pects of this scenario are still controversial, observations 
confirm the connexion between QSO, ULIRGs and merg- 
ers (Canalizo & Stockton 2001, and references therein). In 
conclusion, rather than segregating LIRGs and ULIRGs 
into starbursts and AGNs, the key question is to constrain 
the relative contribution of each component to the IR lu- 
minosity in a case by case basis. 

3.1. Emission-line diagnostic diagrams 

An important constraint on the contribution of AGN ra- 
diation, if any, to the total luminosity of our galaxies comes 
from the use of emission-line diagnostic diagrams. The set 
of emission lines displayed in Table 2 allows us to examine 
two of the most commonly employed diagrams. In Fig. 5 
we compare the ratio [O iii]A5007/H/3 (which is an indica- 
tor of the gas excitation) versus the ratios [N ii]A6583/Hq! 
and [S ii]AA6716,6731/Ha (the tabulated line ratios are 
also given in Table 8). Two interesting properties of these 
diagrams are their ability to properly discriminate between 
mechanisms responsible for the gas ionization (star forma- 
tion, AGN, or shocks), and their low sensitivity to ex- 
tinction corrections. The latter is clearly manifest in our 
data, since the mean correction of the plotted ratios in 
our galaxy sample is only 2%, with a maximum of 7% 
for the [O iii]A5007/H/3 ratio in hd2-264.2. In panels 5(a) 
and 5(b) we display a large collection of emission line ratios 
corresponding to nearby galaxies from Veilleux & Oster- 
brock (1987) and Gallego et al. (1996) — see table caption 
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for details — , wiiereas in panels 5(c) and 5(d) we do the 
same with our galaxy sample. The analysis of this figure 
clearly indicates that the locus spanned by our galaxies is 
perfectly compatible with that exhibited by star forming 
galaxies. 

It is also very interesting to compare the measured line 
fluxes exhibited by our galaxy sample, and in particular 
the ratios [O iii]A5007/H/3 (an indicator of the emitting 
gas excitation), with those of local star- forming galaxies, 
as a function of luminosity. For this purpose, we have 
estimated rest-frame absolute blue magnitudes for each 
galaxy using the available HST y606 and /814 magnitudes 
(note that the HST F606W and F814W bandpasses corre- 
spond approximately to the rcst-framc B band at z ~ 0.4 
and z ~ 0.8, which are precisely quite close to the red- 
shifts exhibited by the present sample). The HST fluxes 
were corrected for redshift, the transmission curves of fil- 
ters F606W and F814W were dc-rcdshiftcd to reproduce 
pscudo B bands for each target, and a residual k-correction 
(ranging from —0.06 to 0.28 mag) was derived for each ob- 
ject, using the averaged spectral energy distributions cor- 
responding to Sab, Sbc and Scd galaxies of Fukugita et 
al. (1995). The resulting magnitudes are listed in Table 9, 
and the excitation diagram is displayed in Fig. 6. For com- 
parison, we have included in the same figure the very well 
characterized local sample of the UCM survey (Gallego 
et al. 1996, 1997), using different symbols to discriminate 
among distinct star forming galaxies: SEN (Starburst Nu- 
clei, originally defined by Balzano (1983), are spiral galax- 
ies that host a nucleus with an important star-forming 
process, with Ha luminosities always higher than 10* L©), 
DANS (Dwarf Amorphous Nucleus Starburst, introduced 
by Salzer et al. (1989), are similar to SEN but at lower 
scale), HHH (H ii galaxies hot spot — see Gallego et al. 
(1996) — , are bright galaxies with a global star-forming 
process, optical spectrum dominated by blue colors, and 
strong emission lines; they show similar Hq luminosities to 
those exhibited by SEN, but with large [O iii]/H/3 ratios), 
DHHH (dwarf H II hot spot, are similar to HHH, but with 
Ha luminosities lower than 5 x 10^ L©), and BCD (blue 
compact dwarf, arc star forming galaxies with the lowest 
luminosity and higher ionization). Finally, Seyfert 2 galax- 
ies (Sy2) are also displayed in Fig. 6 to show the plot region 
covered by AGN-dominated objects. The analysis of this 
figure reveals that the galaxies of our sample are perfectly 
compatible with the locus of local SEN, with the exception 
of hd2-264.2, which based on its higher [O ni]/H/3 ratio, 
is located in the region of local HHH objects. 

3.2. Further evidence 

There is a growing body of evidence indicating that the 
population of IR luminous galaxies detected by the ISO- 
CAM deep surveys are generally dust-obscured starbursts 
(Fadda et al. 2002; Elbaz et al. 2002). Genzel et al. (1998) 
and Laurent et al. (2000) report the presence of a di- 
chotomy in the spectroscopic properties of starbursts and 
galaxies whose luminosity is dominated by gravitational 
accretion: the former show stronger MIR broad emission 
features, with a fast decline in their emission shortward 
of A ~ 6/im, whereas AGNs are characterized by a flatter 
spectrum. Through the analysis of these features, Tran 
et al. (2001) showed that AGN activity is dominant in 



local ULIRGs only for very luminous objects. In particu- 
lar, these authors found that the contribution due to star 
formation to the total far-infrared luminosity is, on av- 
erage, 82%-94% for galaxies with < lO^^-^L© (for a 
cosmology with Hq = 75 km s~^ Mpc~^, and qq = 0.5), 
these numbers decrease to 44%-55% at higher luminosi- 
ties. Considering that all our galaxies are below that 
threshold (see Table 9) our data support the idea that 
star formation is the dominant source of ionization in our 
sample (note that the largest Ljr value of ~ 10^^ Lq cor- 
responds to GSS073_1810, which actually consists of two 
interacting galaxies). 

Since accretion around black holes must lead to the for- 
mation of luminous X-ray sources, specially in the hard 
X-ray domain (less sensitive to extinction affecting ob- 
scured AGNs), the cross-correlation of deep X-ray and 
MIR observations is a direct way to check for the mecha- 
nism responsible for the nebular emission- lines. From the 
deep X-ray survey in the HDF-N performed with Chandra 
(Hornschemeier et al. 2000, 2001; Brandt et al. 2001), we 
find that only one of our four galaxies in this field, hd2- 
264.1, is an X-ray source (see Fig. 5 in Brandt et al. 2001). 
Very recently Fadda et al. (2002) have studied in detail the 
mechanisms that power the X-ray emission of ISOCAM 
galaxies in the HDF-N-I-FF presenting X-ray counterparts. 
They found that hd2-264.1 exhibits an X-ray luminosity 
Lx ~ 2 X lO^^erg s^^ and an X-ray spectral index similar 
to those of the star forming galaxies Arp 220 and Arp 244. 
Although we are not going to repeat all the details here, 
it is important to stress that the latter result alone clearly 
suggests that, in spite of its X-ray detection, a dominant 
AGN contribution in this object can be safely dismissed. 

Further support to the idea that the optical light of the 
galaxies under study is dominanted by hot stars is given 
by the LRIS spectra of the GSS subsample. These have 
enough S/N that a cross-correlation of their spectral en- 
ergy distribution (SED) against well-known SEDs of local 
star forming galaxies is easily done. This task was carried 
out by fitting the continuum regions free of high sky-line 
residuals to the SEDs in the spectrophotometric atlas of 
Kennicutt (1992b) masking out emission line regions. One 
of the best fits is obtained with NGC 4775, a typical Sc 
galaxy. This result, shown in Fig. 7, reveals that typical 
high-order absorption Ealmer lines, a signature of young 
stellar populations, are present in the spectra of these tar- 
gets. 

A final piece of data gives more support to the same 
idea: the measured emission-line widths are not broad. 
The averaged FWHM values (derived from the fits to the 
emission lines with better S /N and without high residuals 
of sky lines) after quadratically subtracting the spectral 
resolution, are (in km/sec): 110 ± 9 (hd2-264.1), 94 ± 9 
(hd2-264.2), 101 ± 7 (hd4-656.1), 129 ± 9 (hd4-795.111), 
216 ± 26 (GSS073_1810A), 167 ± 45 (GSS073.1810B), 
327 ± 42 (GSS084_4515), and 216 ± 35 (GSS084_4521). 
With this result, we are confident that none of the galaxies 
under study is a Seyfert 2. 

In conclusion, although we cannot completely rule out 
AGNs as a possible source of ionizing radiation, all evi- 
dence compiled here suggests that such a contribution, if 
present, is very small. 
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4. MULTIWAVELENGTH STAR FORMATION RATES 

We have transformed the previous flux measurements 
and their corresponding errors (Tables 2. 4 and 5) into 
star formation rates, using well-known calibrations. The 
summary of this computation is presented in Table 7, dis- 
played in Fig. 9, and described in this section. 

4.1. SFR from ultraviolet continua 

Since the resulting UV spectrum of star forming galaxies 
is roughly flat in the range of 1500 2800 A, the SFR can 
be derived from any measurement of the UV luminosity 
in that spectral interval, using a calibration such as that 
given in Eq. (1) of Kennicutt (1998), namely, 

SFRuv(MQ/yr) = 1.4 x IQ-^^ Luv(erg s'^ Hz"!). (1) 

Note that at the mean redshifts of the HDF-N and GSS 
galaxies of our sample, the U band magnitudes displayed 
in Tables 4 and 5 —the HST F300W ((A) ~ 3000 A) and 
the KPNO ?7 ((A) ~ 3670 A) fllters— sample practically 
the same rest-frame wavelength Arest ~ 2050 A. The esti- 
mated luminosities were corrected for internal reddening 
using the recipe of Calzctti ct al. (2000) for starbursts, 
where the color excesses for the stellar continuum were 
derived using E(B-V)cont = 0.44 E(B-V)gas (Calzetti 
1997). In addition, the Galactic extinction has been cor- 
rected using the Eg(B— V) from Table 1 and the extinction 
curve of Fitzpatrick (1999). 

The errors were determined numerically by bootstrap- 
ping the original spectra and their corresponding error 
spectra. Thus, the bootstrap propagates the errors in 
E(B— V)gas quoted in Table 2. However, we have not 
included here the intrinsic scatter in Eq. (1). In fact, 
different calibrations between UV flux and SFR differ by 
~ 0.3 dex (Kennicutt 1998). In addition, systematic errors 
may also be present due to the fact that Eq. (1) is valid 
only for galaxies with continuous star formation over time 
scales of 10^ or longer. 

4.2. SFR from [O ii] and Ha luminosities 

Here we use the average calibrations given in Eq. (3) 
and (2) of Kennicutt (1998), 

SFR[o ii](MQ/yr) = 1.4 x lO'^^ ii](erg s'^), (2) 
and 

SFRHa(M0/yr) = 7.9 x lO'^^ ^^^^(erg g-^), (3) 

with the correction for internal reddening was applied as 
described above. It is important to keep in mind that Lua 
traces the instantaneous SFR, because the Ha emission is 
due to the re-processing of ionizing radiation shortward of 
912 A, which is produced by the most massive stars. Thus 
they are strongly dependent on the adopted initial mass 
function. 

Since the errors in E(B— V)gas and the extinction- 
corrected Ha flux arc correlated, the error spectra were 
also bootstrapped to calculate the errors for SFRho- We 
should also note that in this procedure we have not in- 
cluded the intrinsic scatter in Eqs. (2)- (3). 

Finally the values derived for SFR[o ii] and SFRho , and 
their corresponding error bars, were corrected for aperture 
effects. For that purpose, we used the factors listed in Ta- 
ble 3 for the bluest filters. In particular, HDF-N galaxies 



were corrected using the fraction of galaxy light passing 
through the U^oo filter, whereas for the GSS objects we 
employed the factors associated to the Veoe filter (slit 1.18" 
for SFR[o II]; slit 0.76" for SFRhq:)- These aperture cor- 
rected SFRs are the ones tabulated in Table 7. 

4.3. SFR from mid infrared fluxes 

The SFR of a galaxy may be estimated from its inte- 
grated IR luminosity, Lir = L{8 — 1000 fim) using Eq. (4) 

of Kennicutt (1998), 

SFRiR(M0/yr) = 4.5 x 10-^^ LiR(erg s-^). (4) 

The IR SED of local star forming galaxies generally peaks 
in the FIR, between A ^ 40 and 100 /im. In this spec- 
tral region, most of the luminosity is due to the thermal 
emission of the so-called big dust grains. The MIR flux 
is a combination of broad emission features and hot dust 
continuum as in the SED of the proto-typical starburst 
galaxy M 82 presented in Fig. 8a. Our sample of galaxies 
was observed in the LW3 filter of ISOCAM, centered at 
15 /xm (A = 12-18 /xm). 

At the redshift of the HDF-N galaxies (z ^ 0.4), the 
IS0CAM-LW3 filter (upper thick horizontal segment on 
Fig. 8a) samples almost the; same rest-frame wavelength 
as the IRAS-12/im filter for galaxies at z= 0. In Fig. 8b, 
we show that the 12 /xm luminosity {vL„) of local galax- 
ies is tightly correlated with their integrated IR luminos- 
ity, LiR,. We used the 293 galaxies from the IRAS Bright 
Galaxy Sample (BGS Soifer et al. 1989) with 12 ^um fiux 
densities and redshifts to derive this correlation. 

At the redshift of the GSS galaxies {z ~ 0.8), the 
IS0CAM-LW3 filter (lower thick horizontal segment on 
Fig. 8a) largely overlaps the same rest-frame wavelength as 
the IS0CAM-LW2 filter (upper dashed rectangle), which 
is centered at 6.75 /xm (A = 5-8.5 /xm). In Fig. 8c, we 
have reproduced the Fig. 5d from Elbaz et al. (2002) which 
shows that rest-frame LW2 luminosities are also correlated 
with LiR. 

The previous two correlations suggest that we can com- 
pute LiR for galaxies with z ^ 0.4 and 0.8 which were 
detected at 15 ^m. However, because the rest-frame wave- 
length range corresponding to the observed 15 /im ISO- 
CAM band does not exactly match the IS0CAM-LW2 
and IRAS-12/xm filters, we did not compute ijR directly 
from the correlations but instead from a library of tem- 
plate SEDs ranging from 0.1 to 1000 /im. These templates 
were produced by Chary & Elbaz (2002) who combined 
observed SEDs in order to reproduce the same observed 
trends between mid-infrared and far-infrared luminosities 
derived from ISOCAM, IRAS and SCUBA observations of 
nearby galaxies. 

This computation relies on the hypothesis that the 
trends observed for local galaxies apply also to galaxies 
located at 2; ~ 0.4 and 0.8. As suggested in Elbaz et al. 
(2002), this assumption is supported by the correlation 
observed between radio and MIR luminosities both in the 
local and distant universe (z ~ 0.8). 

To obtain the total integrated IR luminosity for the 
galaxies of our sample, we proceeded as follows. We com- 
puted the expected flux density, S',y(15/xm), for all the tem- 
plate SEDs, at the redshifts of each of our sample galaxies. 
The values thus obtained were compared to the measured 
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ISOCAM-15/im flux densities in order to select, for each 
sample galaxy, the template SED that best matched its 
ISOCAM flux density. The chosen SED was then used to 
predict IRAS flux densities. The latter were employed to 
compute the total IR flux densities, using 

(Sanders & Mirabel 1996), where /12, /25, /eo and /loo are 
the IRAS flux densities in Jy at 12, 25, 60 and 100 /zm. 

Note that we did not compute Lir by directly integrat- 
ing the template SED because the spectral energy distribu- 
tions of Chary & Elbaz (2002) were built to fit correlations 
using Eq. 5. Anyway, the difference in the computation of 
LiR, between using Eq. 5 and integrating the SED is of 
the order of only 5-10%. The final IR luminosity for each 
galaxy is listed in Table 9. 

The error bars in SFRir quoted in Table 7 have been 
computed by propagating the errors in the ISOCAM- 
15 jim flux (i.e. obtaining the best SED which fitted 
5,^(15 /xm)±error). They do not include a possible intrinsic 
scatter within the correlations shown in the Figs. 8b, c since 
the observed scatter combines this intrinsic scatter with 
measurement errors. If the scatter observed in Figs. 8b, c 
were mostly intrinsic, the 1-sigma error bars on SFRir 
quoted in Table 7 should be increased according to a 40% 
uncertainty on Ljr. 

4.4. SFR from, radio fl,uxes 

The correlation between radio and infrared luminosities 
(Helou et al. 1985; Condon et al. 1991; Yun et al. 2001) 
offers an additional diagnostic tool to estimate the SFR. 
This correlation can be characterized with the help of the 
flux density ratio q deflned by Helou et al. (1985), 



log 



F(40-120 /im) 
3.75 X 1012 W m-2 



loe 



5'^ (1.4 GHz) 
W m-2 Hz-i 



(6) 

where 5*1^(1.4 GHz) is the observed 1.4 GHz flux density 
in units of W m~^ Hz~^. In the study of a complete sam- 
ple of 1809 galaxies from the IRAS 2 Jy sample, Yun et 
al. (2001) have shown that 98% of the galaxies follow the 
above relation with q = 2.34 ± 0.01, the scatter in the 
linear relation being ~ 0.26 dex. 

We computed rest-frame radio luminosities at 1.4 GHz 
assuming that oc (as suggested by Yun et al. 
2001), using the radio fluxes given in Table 6, and the 
typical radio spectral index a = 0.80 ± 0.15, with the ex- 
ception of hd2-264.1 (a = 0.87 ± 0.12) and GSS084_4521 
{a = 0.75±0.25), whose values were determined more pre- 
cisely due to the availability of radio measurements at two 
different frequencies. 

To convert the derived FIR luminosity, L(40-120mwm), 
into integrated IR luminosity, L(8-1000mum), we used 
the following tight correlations derived from the IRAS 
BGS (Soifer et al. 1989) 

L(8-1000 Atm) = (1.26 ± 0.15) L(40-500 /im), (7) 

L(40-500 /xm) = (1.50 ± 0.10) L(40-120 /xm), (8) 
and thus 

L(8-1000 Mm) = (1.89 ± 0.26) L(40-120 /xm). (9) 



Finally, the SFRradio is computed from Eq. (4), employing 
the L(8-1000 ^m) luminosity obtained from the combina- 
tion of Eqs. (6) and (9). 

The random errors correspond to the coadded effect of 
random errors in the radio flux and in the spectral index a, 
neglecting the contribution of a possible intrinsic scatter 
in the q parameter. 

5. ADDITIONAL GALAXY PROPERTIES 

Although the galaxy sample examined in this paper is 
small, it is clear from Fig. 1 that it includes a diverse 
collection of objects. In addition to star formation rate es- 
timates, the sample can also be used to derive additional 
galaxy properties, like metallicitics and masses, which will 
provide a more complete description of each object. 

Metallicities can be determined from the measured emis- 
sion line strengths using the empirical line flux ratio tech- 
nique of Pagel et al. (1979), which relates the line ratio 
[O ii]A3727+ [O iii]AA4959,5007 



R. 



23 



H/3 



(10) 



to the oxygen abmidance, 0/H. This calibration, shown 
in Fig. 10a, can be parameterized as a function of the ion- 
ization parameter 

[O III] AA4959. 5007 .... 
= [O n]A3727 ' 
In this work we have adopted the analytic fits of Kob- 
ulnicky et al. (1999) in particular their Eqs. (7)-(9) — 
to the set of photoionization models of McGaugh (1991). 
As seen in Fig. 10a, these calibrations are degenerate, 
there being two widely different oxygen abundances for 
a given pair of R23 and O32 values. Fortunately the di- 
agnostic diagram relating [N ii]A6583/[0 ii]A3727 with 
[N ii]A6583/Ha (see Fig. 10b), has proven to be an ex- 
cellent tool to break this degeneracy, as shown by Contini 
et al. (2002). In particular, galaxies that fall in the upper 
right quadrant of this diagram can be safely assigned to the 
upper branch (high metallicity) of the oxygen calibration 
in Fig. 10a. This is the case for all galaxies in our sam- 
ple, excepting GSS084_4515, which falls in the turnaround 
region of the calibration diagram. Even though the mea- 
sures for this galaxy present large uncertainties, it is clear 
that it exhibits the lowest metallicity. Wc; have also in- 
cluded in Figs. 10a and 10b the UV selected galaxies from 
Contini et al. (2002) to illustrate the reliability of this 
method. The R23 and O32 metallicity indicators, together 
with the final derived metallicities, are listed in Table 8. 
Note that for the two colliding galaxies GSS073_1810 the 
[O lll]AA4959,5007 emission lines fall, unfortunately, out- 
side the observable wavelength range of LRIS, so we can- 
not derive their metallicity. However, the location of these 
two objects in Fig. 10b (shown as open circles), suggests 
that they also should have a metallicity corresponding to 
the upper branch of the R23 calibration in Fig. 10a. This 
result is fully consistent with the position of our galaxies 
in the excitation diagram presented in Fig. 6, in which the 
local starburst trend is typically interpreted as a sequence 
governed by the metallicity content of the ionized gas, with 
metallicity increasing when moving from the H II galax;ies 
to the SBNs. 

The locus of our galaxy sample in the metallicity- 
luminosity diagram is examined in Fig. 11, and compared 
with the regions spanned by additional samples of galaxies. 
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extracted from Fig. 10a of Contini et al. (2002, see refer- 
ences therein). In this sense, hd2-264.1, hd4-656.1, hd4- 
795.111, and GSS084_4521 are perfectly compatible with 
the locus spanned by local SBN galaxies (open circles). In 
addition, GSS084_4515 overlaps with the luminous com- 
pact galaxies (LCGs) of Hammer et al. (2001), whereas 
hd2-264.2 is close to the location of the z ~ 0.1-0.5 
emission-line objects of Kobulnicky & Zaritsky (1999). Fi- 
nally, since the positions of the two interacting galaxies 
GSS073_1810 in Fig. 11 indicate that their metallicities 
are likely to be in the range covered by the remaining spi- 
ral galaxies of the sample, and considering their rest-frame 
absolute B magnitude, their location in Fig. 11 should also 
be compatible with the SBN objects. 

Galaxy masses for the HDF-N subsample have been 
estimated using the K band photometry published by 
Fernandez-Soto et al. (1999) and listed in Table 4. Us- 
ing the spectral energy distribution of late-typc galaxies 
given by Fukugita et al. (1995), we have derived an aver- 
age K-correction of ii'(z) —i^(0) = —0.65 mag. Finally, we 
have employed a mass-to-light ratio 0.93 Mq/L^,©, which 
is the value for the local SBN-like galaxies modeled by Gil 
de Paz et al. (2000). The resulting masses are listed in 
Table 9. 

We also quote in Table 9 values for the galaxy half-light 
radius, which come from the Medium Deep Survey (Rat- 
natunga et al. 1999) fits in the F606W filter for HDF-N 
galaxies (roughly rest-frame B band z ^ 0.4), or from 
GIM2D (Simard et al. 2002) fits to the F814W filter (close 
to the rest-frame B band at ~ 0.8) in the case of GSS 
galaxies. The uncertainties for the HDF-N data were es- 
timated from the variation of the half-light radius values 
when using the F300W and F814W filters. The errors 
for the GSS galaxies combine the uncertainties calculated 
by GIM2D using Monte Carlo simulations, added to the 
variation in the galaxy fits using the F606W images. 

We have combined the computed rest-frame absolute 
B magnitudes, with the half-light radii and the emission- 
line widths, to compare the structural properties of our 
galaxy sample (see Fig. 12) with those of different galaxy 
types (both local and at rcdshifts < 1) compiled and 
pubhshed by Phillips et al. (1997). The unambiguous 
spiral-type galaxies of our sample (hd4-656.1, hd4-795.111, 
GSS073.1810 and GSS084.4521) fall in the region spanned 
by local galaxies of this type in these figures (since these 
objects arc observed face-on, their positions in Fig. 12c 
are only lower limits). The disky galaxy hd2-264.1 is 
also compatible with low luminosity and low mass spi- 
rals, whereas the compact object hd2-264.2 is consistent 
with the CNELG and compact galaxies at higher red- 
shifts of Phillips et al. (1997) and Guzman et al. (1997). 
Interestingly, another compact galaxy from our sample, 
GSS084_4515, behaves in a rather different fashion when 
plotted in these diagrams. The latter result is examined 
in more depth in the next section. 

6. DESCRIPTION OF INDIVIDUAL GALAXIES 

Here we present a summary of the properties of each 
galaxy in our sample. 

6.1. hd2-264.1 



This galaxy appears to be closely connected with hd2- 
264.2 (the difference in radial velocity between both ob- 
jects is only 350 km/s), and the star formation episode 
may well have been triggered as the result of the interac- 
tion between both galaxies. In the HST images hd2-264.1 
shows a central blue condensation and a close to edge-on 
diffuse disk. Based on the standard emission-line ratio di- 
agnostics, this object qualifies as a starburst nucleus of 
moderate luminosity (Mb = —19.7), similar in Ha lumi- 
nosity, far-lR luminosity, color and size to the local galax- 
ies in the same spectroscopic category. 

6.2. hd2-26l2 

This object is a featureless compact galaxy close to hd2- 
264.1. Its appearance in the HST images is highly concen- 
trated, with a half-light radius of ~ 0.6 kpc in the F300W 
band, and ~ 1.5 kpc in the F814W band. Its low lumi- 
nosity and small velocity dispersion indicate a low mass 
system where a star- forming process, likely motivated by 
an interaction with its larger companion hd2-264.1, has in- 
creased its luminosity. Its spectrum confirms those char- 
acteristics typical of a high ionization ongoing star for- 
mation process. The object was not detected by ISO, 
and shows similar properties to the compact star-forming 
galaxies at intermediate redshifts studied by different au- 
thors (Phillips et al. 1997; Guzman et al. 1997; Kobulnicky 
& Zaritsky 1999). 

6.3. hd4-656.1 

This face-on spiral galaxy exhibits two well defined spi- 
ral arms in the HST images. The object shows emission 
lines that correspond to a moderate starburst in the nu- 
cleus. The galaxy properties of this object are similar to 
those of its local counterparts. 

6.4. hd4-795.111 

This galaxy has very similar properties, within errors, 
to those of hd4-656.1. Both show moderate infrared lu- 
minosities and spectral features consistent with a nuclear 
starburst, with a star formation rate of about 2 Mq yr~^ 
as indicated by all the SFR tracers. 

6.5. GSS073_1810 

The HST images of GSS073_1810 reveal a system formed 

by two morphologically disrupted spiral galaxies (Aur 
70 km s~^ in Ha). This large (^ 27 kpc total diameter) 
and luminous system is clearly undergoing a strong burst 
of star formation. It is worth noting that each individual 
component could be considered as being average in size 
and luminosity. The combined system, based on its inte- 
grated IR luminosity, can barely be qualified as an ultralu- 
minous infrared galaxy (ULIRG, log Lir/L0=11.97). The 
ratio SFRiR over SFRna is well above the ~1 observed for 
objects with IR luminosity lower than IO^^Lq, confirming 
a different nature for the object. 

6.6. GSS0844515 

The behavior of this interesting galaxy differs signif- 
icantly from hd2-264.2, the other compact object ana- 
lyzed in this work. It shows the lowest metallicity of the 
galaxy sample, and its value matches that of the sam- 
ple of LCGs at 0.5 < ^ < 0.7 (asterisks in Fig. 11) of 
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Hammer ct al. (2001). Its mass 4 x 10^° Mq, as read 
from Fig. 12c), is also similar to that of the LCG sample 
(10i°-2.5 X 10" Mq). GSS084_4515 has an amorphous 
but compact morphology (half-light radius ~ 2.5 kpc), and 
its spectrum is dominated by a high ionization starburst 
that accounts for the Hll-like emission-line ratios. The 
colors, SFRs, diagnostic diagrams. Ha luminosity, abso- 
lute magnitude and metallicity, all suggest that this galaxy 
is similar to the brightest compact star-forming galaxies 
studied by Guzman et al. (1997), in particular the sub- 
sample at z > 0.7 classified as disk-like. 

6.7. GSSO84.452I 

This object has the appearance of a bright face-on 
spiral galaxy, and its spectroscopic properties are typi- 
cal of systems presenting a nuclear starburst. The ob- 
served discrepancies between SFRir and SFRradio with 
SFRuv; SFRj^Q and SFRho suggest that its nuclear 
starburst is highly obsc;ured. The IR luminosity is very 
high (logLiR = 11.59) but similar to extreme local objects 
such as UCM2250-I-2427. 

7. COMPARISON OF DIFFERENT SFR ESTIMATORS 

Although it can be argued that the calibrations adopted 
in Sect. 4 to derive the SFR estimates listed in Table 7 do 
suffer from large uncertainties, in this work we are pre- 
cisely seeking systematic deviations when using the same 
prescriptions in all the galaxies. The discrepancies found 
by following this approach must then be understood in 
terms of processes that violate the hypothesis of a similar 
star forming scenario in all the galaxies. Among these, 
varying star formation histories and differences in the 
wavelength dependent extinction corrections are expected 
to be two of the most important factors. In addition, the 
detailed modeling of nebular emission from star-forming 
objects by Chariot & Longetti (2001) shows that the sit- 
uation is even more complex since the zero-age effective 
ionization parameter, the gas metallicity and the effective 
dust-to-heavy clement ratio are all parameters which vari- 
ation introduces scatter in the predicted SFRs. 

Even assuming the reduced size of our galaxy sample, 
the graphic comparison presented in Fig. 9 does provide an 
interesting insight concerning how different star formation 
indicators behave in galaxies at intermediate redshifts. 

7.1. SFRijy and SFRua 

The comparison between SFRuv and SFRhq (Fig. 9a), 
shows that for galaxies at 2 ~ 0.4, excluding hd2-264.1, the 
SFRuv values arc slightly above the 1:1 relation. In the 
specific case of hd2-264.1, the HST images suggest that we 
are observing an almost edge-on disk, so the UV extinc- 
tion correction might be underestimated and the actual 
UV flux be slightly higher. Interestingly, all the galax- 
ies at higher rcdshift (z ~ 0.8) exhibit underestimated 
SFRuv- These results are consistent with the work of 
Sullivan et al. (2001), who found that the UV/Ha ratio 
decreases with increasing SFR for galaxies in the redshift 
range < z < 0.5, whereas for low-SFR objects the UV 
luminosities lead to higher SFRs than Ha. Previous ob- 
servations of distant galaxies, using uncorrected Ha fluxes, 
had already revealed that the UV flux underestimates the 



SFR by a factor of a few (Glazcbrook et al. 1999; Yan ct al. 
1999; Moorwood et al. 2000). The origin of this dichotomic 
behavior has tentatively been attributed to temporal vari- 
ations in the star formation histories, with episodic and 
rapid starbursts (Glazebrook ct al. 1999; Sullivan ct al. 
2001). However, Beh & Kennicutt (2001) find this expla- 
nation insufficient, and argue for an increasing attenuation 
of the UV emission relative to Ha for higher luminosity 
galaxies. 

7.2. SFR^Q and SFRna 

Although it is well known that the observed (i.e. uncor- 
rected for extinction) [O Il]/Ha ratio shows a large scat- 
ter in local star-forming galaxies (a factor of ~ 7), this 
scatter improves when extinction corrected values are used 
(Jansen et al. 2001; Aragon-Salamanca ct al. 2002). In our 
case (Fig. 9b) the comparison between both SFR estimates 
is, at first sight, quite good within error bars. However, 
it is interesting to note that the galaxy that deviates the 
most from the 1:1 relation is precisely the object with low- 
est metaUicity, namely GSS084_4515. That metallicity is 
an underlying cause of scatter when comparing [O ii]A3727 
with Ha luminosities has been suggested by Jansen et al. 
(2001). To explore this issue in greater detail, we compare 
in Fig. 13 the correlation found by these authors between 
the [O 11] /Ha ratio and the metallicity parameter R23 for 
local emission-line galaxies which are not dominated by an 
active galactic nucleus. The solid line is the linear fit to 
the local galaxy sample, which is given by (Jansen et al. 
2001) 

log([0 ii]/Ha) = 0.82 log R23 - 0.48. (12) 

The galaxies in our sample follow very well this trend. 
After applying a metallicity correction to the [O 11] SFR 
estimates (Fig. 14) we find an excellent match between 
both SFR indicators. 

7.3. SFRiji, SFRradio and SFRB_a 

Finally, in Fig. 9c and 9d we compare the unbiased 
SFRiR and SFRradio estimators, with SFRho- For the 
two galaxies of our sample with radio detections, SFRir 
and SFRradio agree very well within the error bars. The re- 
maining upper limits in radio fluxes also give consistent re- 
sults with the IR data, excluding GSS073_1810, for which 
the IR flux is unexpectedly higher than the constraint im- 
posed by the radio flux. For the latter galaxy, an inspec- 
tion of the radio map of Fomalont ct al. (1991, Fig. 2, NW 
quadrant) reveals that there is indeed a detection (dis- 
played as an isocontour at a level of 8fj. Jy) at the ex- 
pected location of this object. The fact that it constitutes 
an interacting system with a physical star-forming exten- 
sion which is likely larger than that exhibited by a single 
galaxy like GSS084.4521, might explain why GSS073.1810 
does not present a concentrated peak flux density above 
the 16/x Jy threshold used by Fomalont et al. (1991) to 
generate their "complete" catalog. Thus, the integrated 
flux density for this object is probably very close to the 
sensitivity limit of their survey. This result, together with 
the uncertainty in the total infrared luminosity, lessens the 
discrepancy between SFRir and SFRradio for this object. 

The most intriguing result, which is specially manifest in 
Fig. 9c, is the good agreement between these SFR indica- 
tors for the galaxy subsample at 2; ~ 0.4, while for the 
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higher redshift (and more luminous) objects the agree- 
ment is very poor. To investigate this paradoxical be- 
havior in greater detail we compare in Fig. 15 the ratio 
SFRiR/SFRna as a function of the IR luminosity. In this 
diagram we have included additional galaxy samples, both 
local and at high redshift. The small open symbols cor- 
respond to local star-forming galaxies from Gallego et al. 
(1996) and Buat et al. (2002, excluding cluster galaxies 
and objects with apparent diameters larger than 1.5 ar- 
cmin likely affected by aperture effects ). The small 
solid squares represent the measurements of galaxies ob- 
served by Rigopoulou et al. (2000) in the Hubble Deep 
Field South without a correction for internal extinction. 
By including the average extinction correction of these au- 
thors (a factor of 4 in the Ha flux), their measurements are 
displaced to the positions indicated by the tip of the ar- 
rows. The Ha fluxes in the latter sample were also not cor- 
rected for aperture effects. Wc have ignored two galaxies 
from the Rigopoulou sample, namely ISOHDF-S_38 and 
ISOHDF-S_39, because they are probable AGNs (both ex- 
hibit very strong [N ll] lines compared with Ha — see their 
Fig. 1 — ). Asterisks are the extinction corrected objects 
from Fig. 3a of Sullivan et al. (2001), where we have de- 
termined LiR from their 1.4 GHz luminosities using the 
relations shown in Sect. 4.4 (we have excluded upper lim- 
its in Sullivan et al. data). The dashed line is a bisector 
least-squares fit (Isobe et al. 1990) to our galaxy sample, 

SFRir/SFRh« ~ 3.3 X 10-^ {Li^/ L^f-^" , (13) 

where we have excluded the upper limits in Lir (and thus 
in SFRir/SFRhq:)- Our lower redshift {z ~ 0.4) galax- 
ies are compatible with the local reference sample, in the 
sense that SFRhq traces almost all the SFR (if not all) 
as derived from the IR luminosity. However, this is not 
true when the more luminous galaxies are considered. In 
particular, only one object, out of 10 local galaxies with 
LiR > 10^^ Lq, falls above the relation given in Eq. (13). 
On the contrary, and although the actual extinctions for 
the high-redshift sample of Rigopoulou et al. (2000) are 
unknown, the averaged correction estimated by these au- 
thors seems to accommodate their galaxies in the previ- 
ous fit. Finally, the two most luminous IR objects from 
the sample of Sullivan et al. (2001) do also appear quite 
close to the same fit. Summarizing, the data we have col- 
lected here indicate that the extinction corrected SFRhq 
estimates in z ^ 0.8 galaxies miss an increasing fraction 
of the total SFR, which is clear evidence of the presence 
of highly obscured dust-enshrouded star forming regions 
within these galaxies. In addition, there is an indication 
that this underestimate may be a function of Ljr (Eq. 13). 
However, considering the luminosity and redshift segrega- 
tion of our galaxy sample, we cannot yet preclude an effect 
due only to z. 

Interestingly, the two galaxies with compact aspect, 
hd2-264.2 and GSS084_4515, were both undetected by 
ISO, and thus we have only plotted upper limits for these 
objects in Fig. 15. However, it is important to high- 
light that this common morphology is misleading, since, 
as we have previously discussed, these galaxies exhibit dif- 
ferent metallicity, mass and structural properties, apart 
from a difference of a factor of ~ 10 in total far-infrared 
luminosity. The arrow in hd2-264.2 places this object 
close to SFRiR/Ha ~ 1, whereas in the case of the LCG 



GSS084_4515 its real Lir would put it either on the fit 
given by Eq. (13), or in the locus of the z ^ 0.4 ob- 
jects. If, as suggested by Hammer et al. (2001), LCGs are 
evolved starbursts and the progenitors of the present-day 
spiral bulges, objects of this type may have an important 
contribution to the SFR density at intermediate redshifts. 
These authors have also predicted the total far-infrared lu- 
minosities of LCGs to be around 10^^ L©, near or slightly 
below the limiting sensitivity reached by ISO in the red- 
shift range 0.5 < z < 0.7. If we assume that GSS084_4515 
belongs to the same galaxy family, and although its non- 
detection by ISO is then not surprising, this galaxy will 
accommodate in the fit of Eq. (13). This result highlights 
the relevance of quantifying the actual Lir of LCGs, which 
will be possible in a near future with the help of SIRTF. 

Another important question is whether it is possible to 
use any measurable parameter exclusively derived from the 
optical lines with the aim of correcting the biased SFRrq 
estimates, in order to match the values obtained from the 
total far-infrared luminosity. In this sense, and as we 
already mentioned in Sect. 1, Hopkins et al. (2001) and 
Sullivan et al. (2001) have recently reported a correlation 
between the optical extinction (derived from the Balmer 
decrement) and the star formation rate (or luminosity) of 
the galaxies. Buat et al. (2002) also find correlations be- 
tween the measured extinction in Ha, A[Ha], and Lir, 
although with a very large scatter. 

Wc have represented in Fig. 16 the ratio SFRip /SFRna 
as a function of A [Ha]. For local galaxies this ratio corre- 
lates with A [Ha] when Ha fluxes uncorrected for extinc- 
tion are used. After correcting the Ha flux for extinction, 
local galaxies and our galaxy subsample at z ^ 0.4 scatter 
around SFRir/SFRrq ~ 1, whereas the z ~ 0.8 objects 
remain above this value. Since the measured Ha extinc- 
tions in GSS073_1810 and in GSS084_4521 are very similar, 
wc cannot use this parameter to conclude whether there 
is any correlation between the systematic underestimation 
of SFR from Ha as a function of A[Ha]. Anyway, we 
should note that these two galaxies have Ha extinctions 
similar to hd2-264.1, and lower than hd2-264.2. Both hd2- 
264.1 and hd2-264.2 have total far-infrared luminosities 
one order of magnitude fainter that either GSS073_1810 
and GSS084_4515. For that reason, if there were actually 
a correlation between the observed rest-frame optical ex- 
tinction and the ratio SFRiR/Ha, it should be sought in 
galaxies with ijR > 10^^ Lq and, presumably, at interme- 
diate redshifts. Clearly more data are required in order to 
settle this question. 

8. CONCLUSIONS 

In this work we have compared several star formation 
rate estimators using a small but diverse sample of galaxies 
at two intermediate redshifts {z 0.4 and z 0.8). By se- 
lecting the sample galaxies on the basis of their ISOCAM- 
15 /im micl-infrared luminosities, we have sampled SFRs 
ranging from 2 to 160 M©yr~^. Their morphological types 
include spiral, compact and merging emission-line galax- 
ies. It is important to highlight that the redshift segrega- 
tion of the sample is also accompanied by a segregation 
in total far-infrared luminosity, with the farthest galaxies 
being the most luminous. This fact leads to an ambigu- 
ity between luminosity and distance that must be kept in 
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mind when interpreting the results. 

The observed spectral range includes the most promi- 
nent optical emission lines, from [O ll]A3727 to 
[S ii]AA6716,6731. The availability of and Ha allows 
the determination of color excesses in the nebular gas and 
correcting the emission-line fluxes for extinction. The di- 
agnostic diagrams generated from the corrected emission- 
line ratios confirm the star-forming nature of the galaxies, 
as it was already suggested by evidence from their total 
far-infrared luminosity, observations in the X-ray domain, 
and the absence of broad emission lines. In addition, the 
overall characteristics of the selected galaxies, in particu- 
lar their metallicity and structural properties, also match 
those of star forming galaxies, both local and at interme- 
diate redshifts. 

The good agreement between HST and ground-based 
photometry, in one hand, and the spectrophotometric ca- 
libration of the spectra, on the other, allows us to be con- 
fident on the reliability of the estimated aperture correc- 
tions, the final optical emission-line luminosities, and thus 
the SFRs from them obtained. Our results have shown the 
following. 

(i) There is a general good agreement in the compari- 
son between the five available SFR indicators (from UV, 
[O ll]. Ha, IR and radio luminosities) for the z ^ 0.4 
(LiR < 10^^ L0) galaxies, whereas the situation is not 
the same in the case of the z ~ 0.8 (Lir > 10^-^ Lq) 
subsamplc, being the discrepancies different depending on 
the considered SFR estimator. Focusing in the z ~ 0.8 
galaxies, the SFRs derived from UV, [O 11] and IR lumi- 
nosities are, respectively, lower, similar and higher, than 
the values obtained from Ha. The paradoxical behavior 
of SFRuv versus SFRhq is still subject of debate, and 
among the possible explanations for the observed discrep- 
ancies are temporal variations in the star formation history 
and luminosity dependent attenuations of the UV emis- 
sion. From the observational perspective, we confirm the 
findings of Sullivan et al. (2001) and extend their result 
up to ~ 0.8. 

(ii) The correction for metallicity of the [O il]/Ha ra- 
tio (using the relation found in local galaxies) greatly im- 
proves the concordance between the SFRs derived from the 
extinction corrected luminosities of both emission lines. 
The correction does in fact work for the galaxies in our 
sample in the two redshift (or luminosity) bins. 

(iii) The fit given in Eq. (13), and shown in Fig. 15, 
indicates that extinction corrected SFRho estimates in lu- 
minous infrared galaxies at intermediate redshifts miss an 
increasing fraction of the total SFR, and that the degree 
of underestimation increases with Lir. This result con- 
firms the finding by Rigopoulou et al. (2000) who used an 
averaged extinction correction of ~ 4 derived from V — K 
colors. In this work we quantify more rigorously this effect 
by employing line fluxes corrected for extinction making 
use of color excesses computed from the Balmer decre- 
ment. Thus, there is here an observational evidence that 
the Balmer emission lines, Ha and H/3, do not probe the 
same region of a galaxy than the one which is responsi- 
ble for the strong IR luminosity in the case of luminous 
IR galaxies. This is even more remarquable in Fig. 16 
where galaxies with very different SFRs share the same 
A [Ha]. The spatial resolution in the MIR is not sufficient 



to confirm this statement at these redshifts, but the case 
of the Antennae galaxy NGC 4038/4039 (Mirabel et al. 
1998) offers a perfect example of such behavior in the lo- 
cal universe. A possible explanation for this effect is that a 
fraction of the star formation in these galaxies is embedded 
in dense and opaque dust clouds and that this fraction in- 
creases with the total star formation. As already expected, 
galaxy encounters (hd2-264.1 and hd2-264.2) and mergers 
(GSS073_1810) are the likely causes for triggering star for- 
mation episodes. Whether a quantitative estimation of 
the hidden star formation is definitely inaccessible to rest- 
frame optical spectroscopic studies, and in particular to 
the use of the measured Ha extinction, seems unclear, and 
deserves further research work. 

Some kind of evolution in galaxies with Ljr > 10^-^ L© 
at z < 0.4 cannot be discarded, since, as we have shown in 
Fig. 15, an important fraction of local galaxies with those 
luminosities does not show such a large underestimation 
of the SFR when using Ha fluxes instead of far-infrared 
luminosities. However the scarcity of the galaxy samples 
strongly demands additional observational work in order 
to settle this point. 

The inclusion of two apparently compact galaxies in our 
sample has allowed us to analyze the relation between the 
different SFR estimators applied to this kind of object rel- 
ative to what is measured for spiral-like galaxies. In par- 
ticular, one of these two objects qualifies as a luminous 
compact galaxy, with similar properties to those of the 
LCGs at intermediate redshifts studied by Hammer et al. 
(2001). Surprisingly, this object may be suffering simi- 
lar dust-enshrouded star formation processes as galaxies 
with a very different morphological aspect, like the bona 
fide spiral GSS084_4521 or even the antenae-like colliding 
system GSS073_1810. The confirmation of this still specu- 
lative result, which has important implications concerning 
the accurate estimation of the cosmic SFR density at in- 
termediate redshifts, will be possible after the launch of 
SIRTF. 
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Table 1 
Galaxy sample. 



Galaxy^ 


R.A.^ 


DEC.^ 


^spcc 


EG(B-V)d 


Run<= 


t f 


hd2-264.1 


12 36 49.76 


62 13 13.1 


0.475 


0.012 


4 (0.75") 


1200 


hd2-264.2 


12 36 49.38 


62 13 11.2 


0.477 


0.012 


4 (0.75") 


1200 


hd4-656.1 


12 36 42.91 


62 12 16.3 


0.454 


0.012 


4 (1.25") 


900 


hd4-795.111 


12 36 41.95 


62 12 05.4 


0.433 


0.012 


4 (1.25") 


900 


GSS073_1810 


14 17 42.64 


52 28 45.3 


0.831 


0.007 


1 (1.18"), 5 (0.76") 


3000, 2400 


GSS084_4515 


14 17 40.43 


52 27 19.4 


0.812 


0.008 


3 (1.00"), 5 (0.76") 


3600, 2400 


GSS084_4521 


14 17 40.55 


52 27 13.6 


0.754 


0.008 


2 (1.00"), 5 (0.76") 


3000, 2400 



'^Galaxy identification: for the HDF galaxies we follow Williams et al. (1996), whereas for the GSS 
galaxies we employ a name based on the field and the WFPC chip numbers, and the and first two 
digits of the X and Y centroid coordinates. 

^J2000.0 coordinates for the HDF (drizzled Version 2) and GSS (Rhodes et al. 1994) galaxies. 

^Spectroscopic redshifts obtained from measurements of individual spectral features, and/or cross- 
correlation technique. 

''Galactic color excesses derived from the dust maps of Schlegel et al. (1998) and the software 
available at http://astron.berkeley.edu/davis/dust/. 

^Spectroscopic observing runs with Keck: 1=LRIS (April 1996), 2=LRIS (May 1997), 3=LRIS 
(April 1998), 4=ESI (May 2000) and 5=NIRSPEC (July 2000). Slit widths employed in each case 

arc given within parenthesis. 

'Total exposure times (seconds). 



Table 2 



Measured emission-line fluxes and color excesses. 





Galaxy 


F([0 iij) 


F(H/3) 


771 / \(~\ tttI \ 


F([0 mJ) 


F([N llJ) 


-F(Ha) 


F([N iiJ) 


F{[S iiJ) 


F{[S 11]) 


E(B-V)gas 




AA3727,3729 


A4861 


A4959 


A5007 


A6549 


A6563 


A6583 


A6717 


A6734 




hrl2-2fi4 1 


4.49 (2) 


3.33 (1) 


1 11 CI 'I 


3.99 (1) 


2.46 (1) 


14 7 Cl") 


3.23 (1) 


0.40 (3) 


1.62 (3) 


0.37l°;;i 




0.55 


0.44 


0.72 


0.22 


2.13 


1.4 


0.72 


0.41 


0.98 


hrl2-2fi4 2 


4.98 (2) 


2.72 (1) 


0.21 (3) 




1.27 (3) 


15.0 (1) 


3.61 (1) 




1.00 (3) 


o.56t°;J° 




0.25 


0.25 


0.23 




0.58 


0.7 


0.69 




0.50 


hrl4-fi'ifi 1 


17.9 (2) 


7 4fi n ") 


1 11 C3'l 


3.11 (3) 




22.7 (1) 


6.18 (3) 


4.60 (3) 


2.70 (3) 


n 04+°°® 




0.9 


0.30 


0.61 


1.57 




0.9 


1.60 


1.19 


0.77 


hd4^795.111 


6.29 (2) 


5.58 (1) 


0.35 (3) 


1.21 (1) 


— 


20.4 (1) 


6.59 (1) 


1.98 (3) 


1.57 (3) 


U.^U_g Qg 




0.55 


0.29 


0.16 


0.35 




1.1 


1.34 


0.53 


0.42 


GSS073.1810A 


7.46 (2) 


4.61 (1) 






2.98 (3) 


22.3 (1) 


8.37 (1) 


2.83 (3) 


4.70 (3) 


45+°- 12 




0.10 


0.52 






0.68 


1.3 


1.28 


0.48 


0.78 


GSS073.1810B 


3.87 (2) 


2.92 (1) 






1.65 (3) 


11.2 (1) 


4.88 (1) 


1.73 (3) 


2.03 (1) 


„ 24+0.18 




0.10 


0.48 






0.63 


1.1 


0.90 


0.38 


0.83 


GSS084.4515 


16.4 (2) 


4.35 (1) 


6.18 (1) 


15.3 (1) 


4.20 (3) 


16.2 (1) 


1.83 (3) 


1.00 (3) 


1.20 (3) 


22+° !'' 




0.2 


0.62 


0.76 


0.5 


1.75 


2.2 


1.20 


1.32 


1.55 


GSS084.4521 


5.05 (2) 


3.82 (3) 


0.38 (3) 




2.53 (1) 


17.6 (1) 


10.1 (1) 


5.77 (3) 


4.30 (3) 


„ 41+0.35 




0.09 


1.24 


0.28 




1.88 


1.9 


4.9 


3.08 


2.35 



Note. — Total measured (i.e. uncorrected for intrinsic reddening) emission line fluxes, in units of 10^^'^ erg cm"^. Random errors are 

given below each flux value. Numbers within parenthesis indicate the type of Gaussian fit performed to the emission lines: (1) Gaussian with 3 Q 

free parameters; (2) two Gaussians with the same width, different amplitude, and flxcd separation; (3) Ciaussian with fixed width and position. 9^ 

In the case of the restricted fits (types 2 and 3), the a priory required fixed values have been taken from neighbor and similar lines. [O in] lines D; 

for GSS073_1810 are outside the observed spectral range. The remaining void entries in the table are due to very high sky line residuals. Note £L 

that these numbers have not been corrected for aperture effects. fD_ 
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Table 3 

Aperture corrections estimates. 



Galcixy 


TJ 

slit 0.75" 


-B450 
slit 0.75" 


I'eoe 
slit 0.75" 


slit 0.75" 


hd2-264.1 
hd2-264.2 


0.61 (0.52 ±0.11) 
0.96 (0.75 ±0.30) 


0.58 (0.50 ±0.12) 
0.92 (0.72 ±0.27) 


0.58 (0.50 ±0.12) 
0.90 (0.70 ±0.26) 


0.60 (0.51 ±0.11) 
0.88 (0.66 ±0.27) 




slit 1 25" 


-B450 
slit 1 25" 


^^606 

slit 1 25" 


-^814 

slit 1 25" 


hd4-656.1 
hd4-795.111 


0.85 (0.69 ±0.17) 
0.77 (().(i;-!±0.14) 


0.84 (0.67 ±0.18) 
0.78 (0.()2 ± 0.17) 


0.82 (0.67 ±0.17) 

0.81 (0.6.5 ± 0.17) 


0.81 (0.68 ±0.15) 
0.8() (0.68 ±0.17) 




slit 0.76" 


-^814 

slit 0.76" 


^606 
slit 1.18" 


^^814 
slit 1.18" 


GSS073_1810 


0.52 (0.45 ±0.10) 


0.56 (0.47 ±0.11) 


0.69 (0.65 ±0.06) 


0.72 (0.67 ±0.07) 




^606 

slit 0.76" 


.^814 

slit 0.76" 


^606 

slit 1.00" 


.^814 

slit 1.00" 


GSS084_4515 
GSS084_4521 


0.88 (0.60 ±0.24) 
0.60 (0.50 ±0.10) 


0.90 (0.67 ±0.25) 
0.61 (0.52 ±0.11) 


0.94 (0.76 ±0.22) 
0.73 (0.62 ±0.14) 


0.97 (0.81 ±0.21) 
0.74 (0.65 ±0.12) 



Note. Fraction of light inside the slit for all the observational configurations. Those numbers 
have been measured in the available HST images (F300W, F450W, F606W and F814W for HDF-N 
galaxies, and F606W and F814W for GSS galaxies). The numbers within parenthesis indicate the 
mean and standard deviation in numerical simulations, assuming typical errors in the positioning 
of the slit on the targets of 0.1" in both right ascension and declination, and of 3° in position 
angle. 



Table 4 

Additional flux data for HDF-N galaxies. 



Galaxy ^450^ V^eoe'' /814'' 



hd2-264.1 


24.93 


± 


0.08 


23.52 


± 


0.01 


22.30 


± 


0.01 


21.50 


± 


0.01 


19.93 


± 


0.01 


hd2-264.2 


24.12 


± 


0.03 


23.48 


± 


0.01 


22.70 


± 


0.01 


22.10 


± 


0.01 


21.43 


± 


0.02 


hd4-656.1 


22.84 


± 


0.01 


22.22 


± 


0.01 


21.33 


± 


0.01 


20.74 


± 


0.01 


19.91 


± 


0.01 


hd4-795.111 


23.22 


± 


0.02 


22.55 


± 


0.01 


21.64 


± 


0.01 


21.00 


± 


0.01 


19.92 


± 


0.01 



'^Magnitudes are given in the AB system, where m = —2.5 log — 48.60, with fi, in erg 
cm~^ Hz~^. AH the data extracted from Table 5 of Fernandez-Soto et al. (1999). 



Table 5 

Additional flux data for GSS galaxies. 



Galaxy 


Tj a 


D a 
.^mag 


p a 
-'t'mag 


T 

-'mag 


^606^ 


.^814'^ Jm'^ 


GSS073.1810<= 

GSS084.4515 

GSS084.4.-)21 


23.38 ±0.02 
24..'34 ± 0.04 
24.08 ± 0.02 


23.00 ±0.02 
24.34 ± 0.04 
23.47 ± 0.02 


21.76 ±0.01 
23.24 ± 0.03 

21.85 ± 0.01 


20.98 ±0.01 
22.60 ± 0.04 
21.04 ± 0.01 


22.00 ±0.01 
23.73 ± 0.04 
21.80 ± 0.02 


20.91 ±0.01 20.83 ±0.28 
22.78 ± 0.02 — 
20.98 ±0.01 



^Magnitudes are given in the AB system, where m = — 2.5 log — 48.60, with in erg cm ^ s ^ Hz ^ . Data extracted 
from Brunner et al. (1999). 

''The /mag values have been corrected for the 0.15 mag offset found in the comparison shown in Fig. 4. 

'^HST V606 and 7814 magnitudes in the AB system, measured in this work using SExtractor (Bertin and Arnouts 
1996). 

''NIRSPEC J magnitude (filter N3) measured in this work (see text for details). 
•^Tabulated magnitudes correspond to the integrated values of the two interacting galaxies. 
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Table 6 

mid-ir and radio fluxes for the galaxy sample. 



Galcixy 






hd2-264.1 


115 ±40 


14 ±3 


hd2-264.2 


< 100 


< 9 


hd4-656.1 


49+f 


< 9 


hd4-795.111 




< 9 


GSS073_1810'= 


361 ± 80 


<16 


GSS084_4515 


< 200 


< 16 


GSS084_4521 


240 ± 90 


24 ±4 



'^Mid-infrared flux densities, in /uJy, as 
measured by ISOCAM with the LW3 filter. 

Sec references in the text. 

''Radio flux densities in /xJy, at 8.5 GHz 
for the HDF-N galaxies, and at 5.0 GHz 
for the GSS galaxies. See references in the 
text. 

'^Tabulated fluxes correspond to the inte- 
grated values of the two interacting galax- 
ies. 



Table 7 

SFR ESTIMATES (Mq/yr) CORRECTED FOR EXTINCTION AND FOR APERTURE EFFECTS. 



Galaxy SFRuv A(UV)^ SFRp „] A([0 ii])^ SFRho A(Ha)^ SFRir SFR^adi, 



hd2-264.1 


1 0+0.8 


1 r,+0.6 


4.4l|8 


1 7+0.7 
^- ' -0.6 




0.88t°:^^ 


6.413-3 


11 ±4 


hd2-264.2 


c 1+2.0 
^•-'--1.4 


9 O+0.4 
^•^-0.3 


7 1+3.5 
' -^-2.3 


^•"-0.4 


q R+0.9 
C).D_Q 7 




< 5.2 


< 20 


hd4-656.1 




^•^"-"-O.IS 


9 8 + 0-7 


"•^^-0.23 


1 «+0-2 
-'-•^-0.2 




-^"'-O.S 


< 18 


hd4-795.111 






1 q+0.6 


97+0-31 


9 9+0.4 
^•^-0.4 


"•+^-0.15 




< 16 


GSS073_1810A 






oc+22 
■^"3-13 


9 1+0.5 










GSS073_1810B 








1 1+0.8 
^■^-0.7 




0.58t°1? 






GSS073_1810^ 






42+23 


-r 7+0.3 
l- ' -0.4 




0.8410-2 


1611^5^ 


<86 


GSS084_4515 


Q 1 +2.6 
o-i_1.4 


0.8912;^^ 


r,r,+22 
^"-10 


1 0+°* 
^•^-0.8 


7.4;^-? 


0.53l°:l 


< 67 


<81 


GSS084_4521 


' -0-4.7 




1 c+SO 


1 

i.y_j 2 


15+f 


0.96tojl 


671^^ 


94 ±44 



'^Extinction correction factors (in magnitudes) derived from the color excesses given in Tables 1 and 2, 
and applied to correct the UV, [O ii] and Ha fluxes listed in Table 2 in order to get the SFR values 
quoted in this table. 

''SFRjo ii] and SFRhq values for this object correspond to the sum of the respective values for the 
individual galaxies (listed in the previous two entries). 

Note. — See Fig. 9 for a graphical comparison of these numbers. 



Table 8 



Emission line ratios and metallicities. 



Galaxy 


[O in]A5007 


[N ii]A6583 


[S n]AA6716,6731 


N 11 


A6583 


R23^ 




12 + log(0/H)<= 


H/9 
np 


Ha 


Ha 


O II 


A3727 


hd2-264.1 


1.14 ±0.16 


0.22 ± 0.06 


0.13 ±0.07 


0.32 ±0.10 


3.40 ± 0.55 


0.75 ±0.14 


8.83 ± 0.06 


hd2-264.2 


> 0.21 


0.24 ± 0.05 


0.13 ±0.05 


0.22 ±0.04 


> 3.47 


> 0.09 


8.78 ± 0.06 


hd4-656.1 


0.41 ±0.21 


0.27 ±0.07 


0.32 ± 0.06 


0.31 ±0.08 


3.09 ± 0.29 


0.22 ± 0.09 


8.84 ± 0.04 


hd4-795.111 


0.21 ±0.06 


0.32 ±0.07 


0.17 ±0.03 


0.67 ±0.15 


1.66 ±0.16 


0.20 ± 0.05 


9.00 ±0.02 


GSS073.1810A 




0.37 ±0.06 


0.32 ± 0.04 


0.43 ± 0.07 








GSS073.1810B 




0.44 ±0.09 


0.33 ±0.09 


0.75 ±0.14 








GSS084.4515 


3.41 ± 0.50 


0.11 ±0.08 


0.13 ±0.12 


0.07 ±0.05 


9.48 ± 1.37 


1.03 ± 0.05 


8.33 ±0.11 


GSS084.452ld 


0.29 ±0.23 


0.57 ±0.28 


0.55 ±0.22 


0.83 ± 0.40 


2.34 ±0.81 


0.19 ±0.14 


8.92 ±0.10 



''R23 = ([O ii]A3727±[0 iii]AA4959,5007)/H/3 
''O32 = ([O iii]AA4959,5007)/[O ii]A3727 

'^Errors correspond to the propagation of the emission line fluxes uncertainties when deriving the metaUicitics from the R23 and O32 
indices (see Fig. 10a). However, additional uncertainties associated to the use of this empirical calibration are ±0.15 dex (Kobulnicky 
et al. 1999). 

■^The [O iil]A5007 flux for this galaxy has been determined from [O lli]A4959. 



Note. — All the emission line fluxes have been corrected for extinction. 

O 

£0 
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Fig. 1. — HST images of the whole galaxy sample. Each panel corresponds to 5" X 5". The upper panel of each object is a composite RGB 
image derived from the combination of the flux calibrated images. For the HDF-N galaxies, the average of F300W and F450W was employed 
as B, F606W as G, and F814W as R. For the GSS galaxies only F606W and F814W were available. For that reason, the former was employed 
as B, the latter as R, and the average of both as G. The lower panels are blue-and-red enhanced compositions, in which the reddest and bluest 
regions have been intentionally exaggerated. For that purpose, all the images for each galaxy were scaled to the same total averaged flux per 
unit wavelength prior to the color composition (in this sense, differential color variations are easier to display). Next, images were assigned 
to RGB using the same criterion than with the upper panels, but before combining them, the weighting factors (1, 1/2, 1/4) and (1/4, 1/2, 1) 
were applied to the {R, G, B) set of each pixel which held R > B and R < B, respectively. For the HDF-N galaxies the solid lines indicate the 
orientation and width of the slit in the observation with ESI. For the GSS galaxies the solid and dashed lines correspond to the slit employed 
in the observation with NIRSPEC and LRIS, respectively. We remind the readers that total exposure times (in seconds) for HDF-N images 
are 153700 (F300W), 120600 (F450W), 109050 (F606W), and 123600 (F814W). For GSS073 these numbers are 24400 (F606W) and 25200 
(F814W), whereas the integrations in GSS084 are the shortest, 2800 (F606W), and 4400 (F814W). 
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Table 9 

Additional galaxy properties. 



Galaxy 




Rc (kpc)b 


Mass (M0)'= 




hd2-264.1 


-19.7 


c (:q + 2.63 
O.0y_o 24 


4.8 


X 


IQio 


1010.57 


hd2-264.2 


-19.4 


1 n2+°-*^ 

-■^•"^-0.40 


1.2 


X 


IQio 


< 1010-48 


hd4-656.1 


-20.6 


q OQ + 0.49 

''•''^-0.11 


4.4 


X 


IQio 


1010.13 


hd4-795.111 


-20.1 


9 f;q+0.93 
Z.UO_Q 21 


3.9 


X 


10^° 


IQlOOQ 


GSS073.1810A 


-21.3 


"3- '"-0.08 










GSS073.1810B 


-21.4 


o 00+0.27 
"J-^^-O.lO 










GSS073.1810 


-22.1 










1011.97 


GSS084.4515 


-20.1 


9 4r+0.28 








< 10"-59 


GSS084.4521 


-21.7 


A 24+0.27 








1011.59 



'^Rest-frame absolute magnitudes in the B band (AB system), 
derived as explained in Sect. 5. 

'^Half-light radii (in kpc), computed as explained in Sect. 5. 

'^Galaxy masses (in solar units) derived as explained in Sect. 5 

*^Total IR luminosity, Lir = L {8-1000 i^m) , in solar units. 
These numbers have been derived as explained in Sect. 4.3. 
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Fig. 2. — Flux calibrated spectra for the HDF-N galaxy subsample. Each gray scale intensity identifies a diflferent order in the ESI echelle 
spectra. Very high sky residuals have been masked only for display purposes. The location of the detected emission lines is marked. The 
insets show the fits to each emission line (in these enlargements sky residuals have not been masked). Inside these panels we have also 
plotted the normalized atmospheric transmission (upper curve, with the zero level marked by the short horizontal line in the upper left side), 
the subtracted sky spectrum (lower curve, in arbitrary units), and the resulting emission line fit (thin line overplotted on spectrum). The 
zoomed wavelength range is 1600 km/s (3200 km/s for the Ha & [N ll] panel), with ticks every 10 A (observed wavelength). The broad-band 
photometry collected in Table 4 is also displayed. Error bars in the spectral direction represent the filter coverage, whereas the arrows indicate 
the decrease in fiux after applying the aperture corrections given in Table 3. There is a good agreement in the absolute flux calibration of the 
spectroscopic and broad-band photometry data. The residual fringing present in the ESI data is responsible for the high residuals in the Ha 
and [S ii] regions. 
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Fig. 3. — Flux calibrated spectra for the GSS galaxy subsample. In this case, LRIS spectra are plotted together with NIRSPEC data. See 
caption of Fig. 2 for explanation. The broad-band photometry corresponds to data collected in Table 5. Note that for these objects, the 
agreement in the absolute flux calibration of the spectroscopic and broad-band photometry data is also good, exception for the Jn3 value 
of GSS073_1810. However, in this case the uncertainty is much larger than in the rest of the broad-band measurements (note that we have 
plotted an upper error bar segment in flux for this point). 
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/814 (HST) 

Fig. 4.- Comparison of the KPNO /mag with HST /gi4 magnitudes (AB system). The data correspond to measurements of objects in field 
GSS073 (the GSS field with the largest exposure times — see caption of Fig. 1 — ). KPNO magnitudes come from Brunner et al. (1999). HST 
magnitudes have been measured in this work with SExtractor (Bertin and Arnouts 1996). The three GSS galaxies of our sample are plotted 
as big circles. The dashed line is the 1:1 relation, whereas the solid line is the weighted least-squares fit to a straight line of slope unity. The 
offset between both lines is 0.15 magnitudes. 
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-2 -1 -2-10 

Log([N II]A6584/H«) Log([S II]AA6716,6731/Ha) 

Fig. 5. — Diagnostic diagrams for nearby galaxies — panels (a) and (b) — and for our galaxy sample — panels (c) and (d) — . Panels (a) 
and (b): filled symbols, galaxies with active nuclei; open symbols, galaxies with star formation; small symbols, galaxies from Figs. 1 and 2 of 
Veilleux & Osterbrock (1987, hereafter V087) (we only distinguish between emission due to star formation and due to active nuclei; we have 
also included the objects classified as NELGs by V087 — plotted with asterisks — , which are either LINERs or H II galaxies); large symbols, 
UCM galaxies from Gallego et al. (1996) (complete sample of Ha emission-line galaxies at 2 < 0.045); thick solid line: separation between 
diflferent type of galaxies. There is a small but clear oflFset between the locus of star-forming galaxies in V087's data, when compared with 
the UCM galaxies. This shift is probably caused by the use by V087 of data from the literature, which used different slit apertures. Panels 
(c) and (d): emission line ratios for our galaxy sample after applying extinction corrections (triangles for objects with z ~ 0.4, and circles 
for galaxies with z ~ 0.8); arrows indicate that the [O ni]/H/3 ratio is unknown. In the last two panels we have labeled the symbols with 
the letters C, S, D, and A to indicate the morphological type of each galaxy, as indicated in the symbol key. Since the [O in]A5007 line for 
GSS084_4521 was outside the spectral range of LRIS, we have estimated its value using the measured [O iii]A4959 flux. In the case of the 
colliding galaxies GSS073_1810 both lines were not within the observed spectral interval, so only [N n]A6583/Ha and [S ii]AA6716,6731/Ha 
are indicated. Both panels indicate that all the galaxies of our sample fall in the region of star-forming galaxies. 
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Fig. 6. — Excitation as a function of rcst-framc absolute blue magnitude for the galaxies in our sample (large symbols, labeled as in Fig. 5), 
compared to the sample of nearby UCM (z < 0.045) emission-line galaxies of Gallego et al. (1996, 1997) (small symbols). The rest-frame 
absolute B magnitudes for the UCM sample have been extracted from the data presented by Perez-Gonzalez et al. (2000), and were corrected 
for Galcictic extinction using the dust maps of Schlegel et al. (1998). Different symbols are employed to show the classification of the local 
galaxies of different types (see discussion in Sect. 3). 
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Fig. 7. — Comparison of the GSS subsamplc spectra with a local Sc galaxy spectrum. The thin line (light gray) is the spectrum of each 
galaxy (in arbitrary units); the upper spectrum is the normalized atmospheric transmission; the dotted line is the corresponding sky spectrum 
(in arbitrary units and shifted downwards); the thick line is the spectrum of NGC 4775, a typical Sc galaxic from the spectrophotomctric atlas 
of Kennicutt (1992b), rcdshiftcd and scaled to match the contirmum level of each galaxy spectrum. The wavelength location of the Balmcr 
lines present in the plotted spectral range are also shown. Note that, in spite of the unavoidable large residuals associated to bright sky lines, 
there is good agreement between the GSS spectra and the scaled Sc spectrum. Since in all the panels the SED of NGC 4775 is plotted with 
the same scale, the diflferent strength of [O n]A3727 reveals that the GSS galaxies span a range in the equivalent width of this emission line. 
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Fig. 8. — Panel (a): SED of the proto-typical starburst M 82 (as in Elbaz ct al. 2002). The two thick line segments mark the width 
and position of the ISOCAM-LW3 filter if the galaxy were located at z = 0.4 (upper one) and z = 0.8 (lower one). The dashed rectangular 
regions mark the position and width of the IRAS-12^m filter (upper one) and ISOCAM-LW2 filter (lower one). Panel (b); Integrated IR 
luminosity, Lir, versus IRAS-12/im monochromatic luminosity in uL,^ for 293 IRAS BGS galaxies. Panel (c): Integrated IR luminosity, 
^IR) versus ISOCAM-LW2 (6.75 /im) monochromatic luminosity in uL^. This figure, which contains 91 galaxies, comes from Fig. 5d of Elbaz 
et al. (2002). 
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Fig. 9. — Comparison of the five SFR indicators listed in Table 7. Triangles correspond to the HDF-N galaxies (z --^ 0.4) and circles to the 
GSS subsample {z ~ 0.8). Symbols are labeled as in Fig. 5. Upper limits in SFRij^ and SFRj-adio are shown with arrows. See discussion in 
Sect. 7 
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Fig. 10. — Panel (a): Calibration of oxygen abundance, as a function of the line ratio R23 = {[O n]A3727+[0 ni]AA4959,5007)/H/3, 
and parameterized by the ionization parameter O32 = ([O in] AA4959,5007)/[O n]A3727. The dotted lines correspond to the analytic fits by 
Kobulnicky et al. (1999) to the McGaugli (1991) set of photoionization models. The shaded areas are the galaxies of our sample (lighter regions 
galaxies at z ~ 0.4; for hd2-264.2 we only have lower limits in R23 and O23, which have opposite effects in metallicity, as shown by the arrows). 
The dashed lines mark the turnaround region, where the uncertainties arc large. Data points come from Fig. 4 of Contini et al. (2002): open 
squares and triangles discriminate between the lower and upper branches, whereas the open stars arc objects that fall in the intermediate region. 
Panel (b): Diagnostic diagram employed to break the degeneracy in panel (a). We have also included in this diagram the galaxy sample of 
Contini et al. (2002), which illustrates the method to break the degeneracy of the oxygen abundance as a function of the R23 index: triangles 
correspond to their galaxies that fall on the upper branch of the R23 calibration (log([N 11] A6583/Hq) > —1, and log([N 11] A6583/[0 n]A3727) > 
-1.05), squares to objects that he on the lower branch (log([N ll]A6583/Ha) < -1, and log([N ll]A6583/[0 n]A3727) < -0.8), whereas stars 
correspond to galaxies that fall in the turnaround region of the calibration. As is clear from this figure, the metallicity for the galaxies of our 
sample must be read from the upper branch calibration of panel (a), exception made for GSS084_4515, which lies in the intermediate region. 
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Fig. 11. — Metallicity— luminosity diagram. The different galaxy samples have been extracted from Contini et al. (2002, Fig. 10a). Note that 
in order to perform the comparison, the rest-frame absolute B magnitudes of the galaxies of our paper have been transformed to a cosmology 
with Ho = 100 km Mpc^^, and qo = 0.5. The comparison samples include non-local emission line galaxies from Kobulnicky & Zaritsky 
(1999, 4 CNELGs of their sample are represented with filled stars) and 14 luminous compact galaxies from Hammer et al. (2001). The solid 
line is the least-squares fit to local irregular and spiral galaxies from Kobulnicky & Zaritsky (1999, Fig. 4). 
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Fig. 12. — Comparison between blue luminosities, half-light surfa<;e brightnesses and velocity widths of our galaxy sample and the data 
from Phillips et al. (1997). Panel (a): Rest-frame surface brightness (averaged within the half-light radius) and blue luminosities. The 
dotted regions indicate the plot domain spanned by different types of local galaxies. Panel (b): Comparison of rest-frame half-light radii 
and blue luminosities. Panel (c): Comparison of rest- frame half-light radii and velocity widths. Since the clear spiral-type galaxies of our 
sample are viewed almost face on, we assume that the measured emission-line widths are underestimating the actual rotational velocity of 
these objects. For that reason they arc plotted as lower limits. The dashed lines are iso-mass tracks (in solar units) corresponding to the virial 
mass estimation used by Guzman et al. (1996), M ~ 3c2/G(T^Re, where we have assumed the exponential case C2 = 1.6 for the geometry 
dependent parameter. In this comparison we transformed our data to a cosmology with Hq = 50 km Mpc~^, and go = 0.05. 
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Fig. 14. — Comparison between the SFRs derived from [O ll] and Ha luminosities without any further correction — panel (a) — and 
employing the metallicity correction log([0 n]/Ha) = 0.82 log R23 — 0.48 (Jansen et al. 2001) — panel (b) — shown in Fig. 13. Note that 
although for the antcnae-like object (the pair of galaxies GSS073_1810) we have applied no correction since R23 is unknown, the locus of 
this system in Fig. 10b suggests that its R23 parameter should be close to that of GSS084_4515, hd4-656.1 and hd4-795.111, for which the 
metallicity corrections in this diagram are small. 
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Fig. 15. — Ratio of SFRih, over extinction-corrected SFRna, as a function of Lja — the horizontal scale at the top gives the SFRir as 
computed from Eq. (4) — . Large filled symbols are the galaxies of our sample (triangles for objects with z ~ 0.4, and circles for galaxies with 
z ~ 0.8; symbols labeled as in Fig. 5), with SFRhc< corrected for extinction and aperture as explained in the text. Small filled squares are the 
galaxies observed by Rigopoulou et al. (2000) in the Hubble Deep Field South: the data points correspond to the values without extinction 
correction, while the tips of the arrows indicate the effect of using the average extinction correction (a factor of 4 in the Ha fiux) employed by 
Rigopoulou et al. in their work; the Ha fluxes in this last sample were not corrected either for aperture effects. Asterisks are the extinction 
corrected objects from Fig. 3a of Sullivan et al. (2001), whereas open circles and stars correspond to local galaxies from Gallego et al. (1996, 
1997) and Buat et al. (2002, excluding cluster galaxies and objects with apparent diameter larger than 1.5 arcmin), respectively. The dashed 
line is a bisector least-squares fit to our galaxy sample, excluding upper limits in Lir (in both axis), whereas the dotted lines indicate the 1 <T 
error in the fit prediction derived from numerical simulations via error bootstrapping. 
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Fig. 16. — Ratio of SFRir over SFRna (with Ha not corrected for extinction in the left panel, and corrected in the right panel) as a 
function of the measured extinction in Ha. Symbols as in Fig. 15. 



